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We present a simple, closed form expression for the potential of an axisymmetric 
disk of stars interacting through gravitational potentials of the form V{r) = —(3/r + 
7r/2, the potential associated with fundamental sources in the conformal invariant fourth 
order theory of gravity which has recently been advanced by Mannheim and Kazanas as 



^ ! a candidate alternative to the standard second order Einstein theory. Using the model 

H : 

. 5t , we obtain a reasonable fit to some representative galactic rotation curve data without the 
need for any non-luminous or dark matter. Our study suggests that the observed fiatness 
of rotation curves might only be an intermediate phenomenon rather than an asymptotic 
one. 
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(I) Introduction 

In formulating a physical theory it is necessary to both work up from phenomenological 
observations and down from fundamental principles, and to be prepared to revise the 
insights obtained from both approaches as new data come on line. However, after a program 
such as this has been successfully carried through once, there is then some reluctance on 
the part of the community to have to reopen the issue even in the light of subsequent data. 
Consequently, the prevailing view on galactic rotation curve data is that their deviation 
from the behavior expected on the basis of the standard Newton-Einstein theory as applied 
to the observed galactic luminous matter surface brightness distribution must be attributed 
to a (rather substantial) non-luminous or dark matter galactic component. Since there is no 
clear evidence today that the dominant component of the Universe is in fact non-luminous, 
there is thus some merit in going back over familiar ground to see where, if anywhere, 
something could be modified. (Milgrom (1983) and Sanders (1990) have also looked at 
possible revisions to the standard theory, with the recent study of Begeman et al (1991) 
in particular showing that Milgrom's MOND alternative is currently phenomenologically 
viable). 

Noting that there is currently no known theoretical reason which would select out 
the standard second order Einstein theory from amongst the infinite class of (all order) 
covariant, metric based theories of gravity that one could in principle at least consider, 
Mannheim and Kazanas have reopened the question of what the correct covariant theory 
of gravity might be (Mannheim and Kazanas (1989, 1991, 1992), Mannheim (1990, 1992, 
1993), Kazanas and Mannheim (1991)), and developed an approach which works down 
from an additional fundamental principle, namely that of local scale or conformal invari- 
ance, the invariance now believed to be possessed by the other three fundamental strong, 
electromagnetic and weak interactions. This invariance forces gravity to be described 
uniquely by the fourth order action Iw = — ck / d'^x{—gY^'^C\^y^C'^'^^'^ where Cx^^^ is 
the conformal Weyl tensor and a is a purely dimensionless coefficient. In their original pa- 
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per Mannheim and Kazanas (1989) obtained the exact, non-perturbative exterior vacuum 
solution associated with a static, spherically symmetric gravitational source such as a star 
in this theory, viz. 

-^00 = l/grr = l-(5{2- 2,(5^) /r - 2,(3^ + 7r - fcr^ (1) 

where {3, 7, and k are three appropriate dimensionful integration constants. As they noted, 
this solution contains the familiar exterior Schwarzschild solution and thereby yields the 
standard exterior Newtonian potential term and its standard general relativistic correc- 
tions whenever the two additional potential terms in Eq. (1) may be ignored. The theory 
thus contains the same solution as the standard Einstein theory in the appropriate kine- 
matic regime even while not containing the Einstein equations themselves, this being all 
that observation can require, and thus nicely meets the demands of solar system scale 
observations. The quadratic term in Eq. (1) may be associated with a general cosmo- 
logical background de Sitter geometry and is otherwise uneventful, and thus (with both 
the 3/?7 terms being numerically negligible - see below) the conformal theory leads to the 
non-relativistic gravitational potential V{r) = —f3/r + 7r/2, which may then be fitted to 
data whenever the weak gravity limit is applicable. V{r) is thus the potential obtained in 
coming down from a fundamental principles approach. In this paper we shall study the 
implications of this potential by working up from galactic data, the first distance scale on 
which the potential is found to lead to significant deviations from Newton. As we shall 
see, the two approaches even have a chance to converge; however, those readers who may 
not be too comfortable with (or even disapprove of) the whole general conformal gravity 
program can view this paper purely as an attempt to identify which phenomenological 
potentials the currently available observational data actually permit. 
(2) The Potential of an Extended Disk 

For an explicit application of conformal gravity to objects such as galaxies, the treat- 
ment completely parallels the Newton-Einstein treatment, i.e. the potential of each of 
the fundamental stellar constituents of the galaxy is first found by solving the relativis- 

3 



tic theory exactly, with the potential of the galaxy then being found perturbatively by 
adding up the potentials of these stars in the weak gravity limit. Thus for our purposes 
here the requisite stellar potential is the above V{r) = — /3/r + with 7 being a new 
stellar parameter about which almost nothing is currently known (7 is however known to 
in principle be completely independent of (3 (Mannheim and Kazanas (1992))), with the 
phenomenological fitting to be presented below providing some first insight into typical 
numerical values for 7. 

In order to actually determine the weak gravity potential of an extended object such 
as a disk, we have found it convenient to generalize the approach of Toomre (1963) first 
to non-thin Newtonian disks (a step also taken by Casertano (1983)), and then to disks 
with linear potentials. The method which leans heavily on the completeness properties 
of the Jn{kR) Bessel functions will be reported elsewhere, and here we state only the 
relevant results. The Newtonian potential of a general axisymmetric distribution of stellar 
matter sources each with the same average /? (i.e. ignoring any morphological variation 
throughout the galaxy and treating each star as though it has the same typical /?) and 
with luminosity density function p{R, z') is calculated to take the form 

nOO nOO nOO 

Vf3{r,z) = -2n(3 dk dR dz' Rp{R, z')jQ{kr)jQ{kR)e-^\^-^'\ (2) 

Jo Jo J -00 

where i?, z' are cylindrical source coordinates and r and z are the only observation point 
coordinates of relevance. Analogously, the net contribution of a disk of stars each with a 
linear 7r/2 potential (i.e. ignoring any possible morphological variation in 7) is found to 
take the simple form 

V^{r,z) = TT^ J dkdRdz' Rp{R,z') 

X [(r^ + R^ + (z- z'f)Jo{kr)Jo{kR) - 2rRJi{kr)Ji{kR)]e~^\^-'''\ (3) 

For a thin exponential disk with p{R,z') = T,{R)5{z') = T,oex.p{—aR)5{z') where Rq = 
1 1 a is the disk scale length and N = 27rEo-Ro ^^e total number of stars in the disk, it is 
possible to perform the integrations completely to yield 

Vf3{r) = -7r/?Eor[/o(ar/2)Ki(ar/2) - h{ar/2)Ko{ar/2)] (4) 
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for the Newtonian contributfon (Freeman (1970)), and 



V^{r) = 7r7Eo{(r/a2)[/o(«r/2)Ki(ar/2) - h{ar/2)Ko{ar/2)] 



+ {r^/2a)[Io{ar/2)Ko{ar/2) + h{ar/2)Ki{ar/2)]} 



(5) 



for the net hnear contribution. Differentiation then yields the remarkably compact expres- 
sion 



an expression which behaves asymptotically as NP/r + N^r/2 — SN'jRl/Ar as would be 
expected. The coefficient NP is usually identified as MG/c^ with M being taken to be the 
mass of the disk; and we see that for thin disks all departures from the standard Freeman 
result are embodied in the 7-dependent term in the simple and compact manner indicated. 
(3) Exponential Disks and Flat Rotation Curves 

In a recent comprehensive analysis of the HI rotation curves of spiral galaxies (the 
more prominent HII optical data studies of the type pioneered by Rubin et. al. (1978) do 
not go out to a large enough number of scale lengths to show any substantive deviation from 
a standard luminous Newtonian behaviour (Kalnajs (1983), Kent (1986)) making the HI 
data the main probe of the outer reaches of the rotation curve), Casertano and van Gorkom 
(1991) have found that the data fall into essentially four general groups characterized by 
specific correlations between the maximum rotation velocity and luminosity; with the four 
groups being intermediate, compact bright, large bright, and dwarf galaxies. Thus as 
a first attempt at data fitting we have chosen to study one representative galaxy from 
each group, specifically the galaxies NGC3198, NGC2903, NGC5907, and DD0154. This 
will immediately enable us to test the flexibility of our theory, as well as confront the 
systematics apparent in dark matter fits to the same four groups where it is typically 
found that the more luminous the galaxy the proportionately less dark matter seems to 
be needed. 



rV'{r) = {N/3a^r^/2)[Io{ar/2)Ko{ar/2) - Ii{ar/2)Ki{ar/2)] 



+ {N^ar^/2)h{ar/2)Ki{ar/2) 



(6) 
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The intermediate galaxy NGC3198 is particularly well-suited for testing theories since 
for it the data go out to the largest known number of surface brightness scale lengths; 
and, with the data being so flat, this galaxy is generally regarded as being prototypical. 
To model the galaxy we have followed van Albada et. al. (1985), and represented the 
surface brightness by a single exponential with a 1' (=2.72 kpc) scale length. (This choice 
approximates Wevers et. al. (1986) U', J, and F filter data with eyeball slopes of Rq = 
63", 58", 54" respectively at a 5% uncertainty level (the F filter data have also been 
confirmed by Kent (1987)), while ignoring a spike in the very small angle region data, and 
also a possible truncation at the edge of the visible region). Moreover, the model ignores 
any modifications to the luminosity profile due to extinction or galactic dust infrared 
reprocessing. Following Begeman (1987, 1989) we have assigned a 2;— thickness to the 
disk according to the general analysis of van der Kruit and Searle (1981), so that the 
luminosity density function p{R, z') needed for Eqs. (2) and (3) takes the separable form 
'E{R)sech.'^{z' / zo)/2zo with zq = Ro/5. (The 2;— thickness structure of the disk is only 
significant at small radii where it serves to ensure that the inner part of the rotation 
curve is well fitted by the Newtonian contribution, thus making it possible to explore 
fully the effect of the linear potential on the outer region). Recognizing a 15% or so 
contribution to the visible mass density from the HI gas itself, we have also included the 
gas as a matter source, and have found that, for model purposes, Wevers et. al. (1986) 
HI surface density data can be well represented by a sum of three exponentials, viz. 
(THiiR) = (37.0 exp(-i2/2.23) + 34.6 exp(-i?/0.87)- 68.2 exp(-i?/1.21)) Mg/pc^ (Ris in 
arc minutes) with a total HI mass (to infinity) of 5.2 x 10^ Mq, of which 4.9 x 10^ Mq 
is observed in the explored 12' region. Finally, again following Begeman, and also van 
Albada and Sancisi (1986), we have multiplied the HI gas profile by a factor of 1.4 to 
allow for the presence of helium. With the model thus defined with only the two free 
parameters and 7 of the stars, we have generated the fit of Fig. (1) to Begeman's 
(1989) rotation curve data. (The fitting proved not to be sensitive to any thickness for the 
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gas so we did not use one, while the fit was also insensitive to any deviations of the value 
of the gaseous 7//? ratio from the stellar one). With N essentially being constrained by 
the overall normalization of the stellar contribution, our best fit is found to yield a value of 
3.8 X 10-^° Mq for the mass of the stars which is quite reasonable for a galaxy with quoted 
luminosity Lb = 9.0x10^ Lbq (i.e. galactic mass to light ratio M/L = 4.2 Mq/Lbq), with 
the obtained galactic mass in fact being a typical so called maximum disk mass in which 
the Newtonian term gets to be as large as it possibly can be. Additionally, we find that 
the coefficient of the net galactic linear term is given by 1/ ^galaxy = ^/^Istar = 2.9 x 10^^ 
cm to yield a galactic gamma to light ratio '^galaxy/ Lb = 3.9 x 10"^*^ /cm /Lbq- As 
we can see from Fig. (1), the contribution of the linear potential piece is remarkably 
reminiscent in shape to that of a typical dark matter contribution to galactic data fitting 
(see e.g. Kent (1987) for an extensive study); and intriguingly we find that the linear 
potential is competitive with the Newtonian one in a galaxy when 1/jgaiaxy is of order the 
Hubble radius, the naive anticipation of Mannheim and Kazanas (1989) in their original 
study. Given this value for ^galaxy, the inferred value for ^star is then 0.9 x 10~^° cm~-^, 
making the linear potential indeed negligible on solar system distance scales as initially 
required, with the linear potential only first becoming competitive with the Newtonian 
one galactically. (In passing we note that with such a small value for 7 the /?7 product 
terms in Eq. (1) are then rendered completely insignificant, a fact we had indicated 
earlier). As regards an assessment of the quality of our fitting, we should note that there 
are some still not fully understood discrepancies in the outer region (of order up to 7 
km/s) between Begeman's data and Bosma's earlier 1978, 1981 data; while additionally 
Begeman quotes a maximum difference of 3 km/s between his inner rotation curve and 
that of Hunter et al (1986). Moreover, we should also note that Begeman's last 2 data 
points (the farthest) actually use adopted values which are extrapolated from closer in 
ones. Also, there is even some indication in the data of a warp at the largest observed 
distances which we have not attempted to model. Finally, in general as regards rotation 
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curve data, we note that velocities are measured at corresponding distances on the two 
sides (receding and approaching) of the galaxy, with differences between the obtained 
values usually constituting the (sole) quoted errors, even though such differences could 
well indicate that the galaxy actually has a symmetry lower than that of a disk (in which 
case the very extraction of the original mass density itself from the surface brightness data 
becomes slightly suspect). Thus fitting at the 5% (or even 10%) level would seem to be 
acceptable, making our fitting quite adequate. 

For the compact, bright galaxy NGC2903 Wevers et al (1986) find a stellar disk scale 
length of 2.0 kpc (=1.08'). For this galaxy the contribution of the gas is quite small 
{Mhi = 2 X 10^ Mq), and introducing a stellar or a gaseous thickness was found to have 
no appreciable effect. Consequently we are able to fit the galaxy directly with Eq. (6) 
(using the same 7//? ratio for stars and gas) to find the fit of Fig. (1) to Begeman's 
(1987) rotation curve data data. We find a fitted mass of 5.3 x 10-*^° M© for the stars 
to be compared with quoted luminosity Lb — 1-5 x 10^^ Lbq {M/L = 3.5 Mq/Lbq), 
and a net galactic linear term given by 1/jgaiaxy = ^/N^star = 1-3 x 10^^ cm, making 
Istar = 1-4 X 10~^° cm~-^ and '^galaxy/ Lb = 5.1 X 10~^° /cm /Lbq- As regards the fitting 
we note that this galaxy also has a warp at large distances which we have not attempted 
to model. Also, in his original dark matter fit Begeman (1987) found a lot of scatter in 
the inner rotation curve (he actually settled for an eye-ball fit rather than a least squares 
one) , prompting him to suggest that there might be an additional stellar component in the 
not well explored close in region. If there is, then in the present theory such a component 
would also contribute in the outer region because of the linear term. 

Until very recently the large, bright galaxy NGC5907 had actually had a severe fitting 
problem because the original surface brightness data of van der Kruit and Searle (1981) 
were simply completely incompatible with the rotation curve data (van Albada and Sancisi 
(1986)) in the inner region where the stellar Newtonian potential should dominate. This 
situation has only recently been rectified by Barnaby and Thronson (1992a) who find (using 
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a different filter) a completely different surface brightness curve, one which does nicely fit 
the inner region (Barnaby and Thronson (1992b)). Barnaby and Thronson find that the 
stellar disk is well parametrized by p{R, z') — T,oexp{—R/Ro)sech.{z'/zo)/nzo where Rq 
=4.0 kpc (=1.22'), and Rq/zq = 9.2. Additionally they find a close in (and thus easy 
to miss) central region stellar bulge with a luminosity 5% of that of the disk which they 
parametrize by a modified Hubble profile (~ l/(-R^ + Rq)) with scale length Rq =0.07'. 
The total stellar contribution to the rotation curve is exhibited in Fig. (1). We find the 
mass of the disk to be 1.1 x 10^^ Mq and that of the bulge to be 1.7 x 10^° Mq with 
the quoted luminosity of the galaxy being Lb — 1.8 x 10^^ Lbq {M/L = 6.1 Mq/Lbq 
for the disk). For the linear term we find (we use the same 7//? ratio for the disk and 
bulge) I /^galaxy = '^/N^star = 1-7 X 10^^ Cm, making ^star = 5.5 X 10~^^ cm-^ and 
Igaiaxy/Ls = 3.2 X 10~^° /cm /LsQ, to givc values which are comparable with those of 
the other galaxies. 

For the dwarf irregular DD0154 Carignan and Freeman (1988) and Carignan and 
Beaulieu (1989) have determined both the rotation curve and the surface brightness data. 
The stellar component is fit by a disk with scale length 0.43' (corresponding to Rq = 0.5 
kpc if the galaxy is at a distance D=4 Mpc - it may be at D=10 Mpc, see below); while 
the gas is well fitted by aniiR) = (31.6exp(-i?/1.42) - 25.7exp(-i?/1.08)) Mo/pc^ {R 
is in arc minutes) with a total HI mass (to infinity) of 2.8 x 10^ Mq (at D=A Mpc), 
93% of which is observed in the explored region. Since the gas turns out to be the main 
gravitational component we see that the observed region corresponds to about 4.5 leading 
1.4' gas scale lengths. In the central galactic region the stellar surface brightness is not yet 
fully explored (it actually appears to be fiattening off, meaning that an exponential disk 
could be overestimating the inner surface brightness) . Consequently in order to fit the inner 
rotation curve we have found it necessary to give the stellar component a secli^{z' /zq)/2zq 
thickness with zq = i?o(stellar)/3. Additionally, we have allowed the stellar and gaseous 
7//3 ratios to vary independently (there is no immediate reason why they should be the 
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same, anymore than the mass to hght ratios of bulges and disks of a given galaxy should 
be the same). We have not included the contribution of gas pressure (a 5% or so effect 
which is completely ignorable in galaxies where stars dominate the dynamics) , nor have we 
considered any random motions of the gas (a 1 or 2% effect), or an apparent galactic warp. 
Also we note that the very last rotation curve data point was only determined on one side 
(the receding one) of the galaxy and its inferred velocity is not independent of those of the 
two immediately previous points. Our best fit at D=4 Mpc is presented in Fig. (1) with 
stellar mass 6.8 x 10^ Mq for a quoted luminosity at that distance of Lb = 5.0 x lO^Lg© 
(M/L = 1.4 Mq/Lbq). For the stars we find l/-f galaxy = l/N-fstar = 4.0 x 10^9 cm, 
making "jstar = 3.7 x 10~^^ cm~^ and ^galaxy/ Lb = 5 x 10~^^ /cm /Lbq, while the best 
value for jgas is found to be zero (i.e. much smaller than the stellar contribution). Now 
while Carignan and coworkers favor putting DD0154 at D=A Mpc, we note that they also 
indicated a possible adopted distance at D=10 Mpc which Krumm and Burstein (1984) 
favor, this being the distance at which the TuUy-Fisher relation is obeyed. Since the amount 
of gas is 6.25 times bigger at the larger distance, it would then be totally dominant. (While 
we shall continue to use the same (jhi{R) as before after scaling up to the larger distance, 
we note that Krumm and Burstein actually obtained a leading scale length of 2.5'). At 
D=10 Mpc the fit is found to be insensitive to any stellar or gaseous disk thickness, with 
the thin disk approximation then yielding the fit of Fig. (1). We find the stellar mass to 
be 7.4 X 10^ -^0- Additionally we find that for the stars 1/ ^galaxy — '^/^Istar — 

1.4 X lO'^" 

cm, making ^star = 9.9 x 10~^^ cm~^, while the best value for ^gaa is again found to be 
zero. (Actually, the fits at both the candidate distances can even be improved in the outer 
region if we allow 7gas to go negative. This would be somewhat difficult to understand, 
though without a dynamical theory for 7 it cannot yet be excluded) . Since our theory, the 
standard ffat dark matter theory, and MOND would all eventually overshoot the data if 
the suggested large distance decline in the DD0154 rotation curve were to be confirmed, 
further observational study of this point might prove interesting. 
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In our fitting we see a reflection of the general luminosity trend found in dark matter 
fits, with both the relative Newtonian contribution and the inferred galactic mass to light 
ratios increasing with luminosity. Additionally, and intriguingly we find that the values ob- 
tained for l/^gaiaxy ^re remarkably close to each other and to the Hubble radius suggesting 
some possible common underlying dynamics. (Perhaps j galaxy sets or is set by the scale 
at which galaxies can fluctuate out of an initial cosmological background). For the three 
regular galaxies NGC3198, NGC2903, and NGC5907 we flnd that the mass to light ratios 
are comparable to each other and likewise their gamma to light ratios, to thus suggest only 
a mild morphological dependence to the average stellar (3 and 7 parameters used as input 
for Eqs. (2) and (3). (While the mass to light ratio is assumed to be uniform within a given 
galaxy, the actual numerical value of this ratio is not apparently universal for all galaxies 
suggesting that galaxies do not all have the same universal mix of stars and/or the same 
typical average stellar p. A similar situation should thus be expected to obtain for the 
galactic gamma to light ratio). We note that the irregular galaxy DD0154 does not show 
the same galactic gamma to light ratio as the other three galaxies. It is not clear whether 
this is a fundamental issue or whether perhaps the galaxy has an anomalous mix of stars. 
As regards this issue we recall that both the dark matter and MOND flts of Begeman et al 
(1991) to the same galaxy flnd atypically small mass to light ratios. Thus dwarf irregulars 
may be fundamentally different, though of course trying to extract out stellar parameters 
in a galaxy where the stars do not dominate may not be completely reliable. With regard 
to our fltting, we see that while our fltting is yielding flat rotation curves, it is doing so 
in a theory in which rotation curves must eventually rise. This stands in marked con- 
trast to the asymptotically flat behaviour expected both in MOND and in the isothermal 
gas sphere model of dark matter with its asymptotically logarithmic galactic potential. 
(Other, less popular, asymptotic alternatives for dark matter are considered in van Albada 
et al (1985)). Since one may unfortunately run out of galaxy before possibly seeing any 
such rise, perhaps the sharpest difference between linearly rising and logarithmic potentials 
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may emerge at the slightly larger distance scale associated with gravitational lensing where 
such differences might even be pronounced, thus making a study of the (so far unknown) 
conformal gravity lensing predictions in the non- asymptotically flat geometry of Eq. (1) 
quite urgent. 

We believe that we have thus established the candidacy (at least) of fourth order 
gravity by working up from the non-relativistic limit; and since the conformal theory has 
already been shown to possess no flatness problem (Mannheim (1992)) and thus not require 
any cosmological dark matter, and since the linear potential has also been shown to be 
capable of yielding galactic stability without the need for dark matter (Christodoulou 
(1991)), we see that both cosmologically and galactically it might turn out to be the case 
that luminous matter is the major constituent of the Universe after all. 

The author would like to thank D. Kazanas, D. Christodoulou and J. Taylor for 
stimulating discussions, and D. Barnaby, K. Begeman, C. Carignan, S. Casertano, S. Kent 
and R. Sancisi for some very helpful communications. This work has been supported in 
part by the Department of Energy under grant No. DE-FG02-92ER40716.00. 
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Figure Caption 

Figure (1). The calculated rotational velocity curves associated with the metric of Eq. (1) 
for the four representative galaxies, the intermediate sized NGC3198, the compact bright 
NGC2903, the large bright NGC5907, and the dwarf irregular DD0154 (at two possible 
adopted distances). In each graph the bars show the data points with their quoted errors, 
the full curve shows the overall theoretical velocity prediction (in km/s) as a function of 
distance (in arc minutes) from the center of each galaxy, while the two indicated dotted 
curves show the rotation curves that the separate Newtonian and linear potentials of Eq. 
(1) would produce when integrated over the luminous matter distribution of each galaxy. 
No dark matter is assumed. 
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